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Introduction to Matrix Algebra 


Definitions: 

A matrix is a collection of numbers ordered by rows and columns. It is customary 
to enclose the elements of a matrix in parentheses, brackets, or braces. For example, the 
following is a matrix: 

5 8 2\ 

-i o V 

This matrix has two rows and three columns, so it is referred to as a "2 by 3" matrix. The 
elements of a matrix are numbered in the following way: 
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That is, the first subscript in a matrix refers to the row and the second subscript refers to 
the column. It is important to remember this convention when matrix algebra is 
performed. 

A vector is a special type of matrix that has only one row (called a row vector) or 
one column (called a column vector). Below, a is a column vector while bis a row 
vector. 


a = 


fy\ 

2 

v3y 


b = (-2 7 4) 


A scalar is a matrix with only one row and one column. It is customary to denote 
scalars by italicized, lower case letters (e.g., x), to denote vectors by bold, lower case letters 
(e.g., x), and to denote matrices with more than one row and one column by bold, upper 
case letters (e.g., X). 

A square matrix has as many rows as it has columns. Matrix A is square but 
matrix B is not square: 


A = 


A symmetric matrix is a square matrix in which x { - = x 7 for all i and j. Matrix A is 
symmetric; matrix B is not symmetric. 
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A diagonal matrix is a symmetric matrix where all the off diagonal elements are 
0. Matrix A is diagonal. 
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9 0 
0 6 
0 0 


0) 
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An identity matrix is a diagonal matrix with Is and only Is on the diagonal. The 


identity matrix is almost always denoted as I. 

'l 0 
0 1 
0 0 


0) 
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Matrix Addition and Subtraction: 

To add two matrices, they both must have the same number of rows and they both 
must have the same number of columns. The elements of the two matrices are simply 
added together, element by element, to produce the results. That is, for R = A+ B, then 


a .i +i ; 
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for all i and j. Thus, 

9 

6 0, 

Matrix subtraction works in the same way, except that elements are subtracted instead of 
added. 


f l 
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Matrix Multiplication: 

There are several rules for matrix multiplication. The first concerns the 
multiplication between a matrix and a scalar. Here, each element in the product matrix is 
simply the scalar multiplied by the element in the matrix. That is, for R = aB, then 


for all i and j. Thus, 


8 


r 

v 


65 

7, 


ab : , 




16 

24 


48A 

56, 


Matrix multiplication involving a scalar is commutative. That is, a B = B a. 

The next rule involves the multiplication of a row vector by a column vector. To 
perform this, the row vector must have as many columns as the column vector has rows. 
For example, 

( 


(1 7 5) 


vC 


is legal. However 
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(1 7 5) 


is not legal because the row vector has three columns while the column vector has four 
rows. The product of a row vector multiplied by a column vector will be a scalar. This 
scalar is simply the sum of the first row vector element multiplied by the first column 
vector element plus the second row vector element multiplied by the second column 
vector element plus the product of the third elements, etc. In algebra, if r = ab, then 


Thus, 




(2 6 3) 1 =2*8 + 6*1 + 3*4 = 34 


All other types of matrix multiplication involve the multiplication of a row vector 
and a column vector. Specifically, in the expression R = AB, 

r = a b 

i) 

where a,-. is the z'th row vector in matrix A and b., is the jth column vector in matrix B. 
Thus, if 

« _in o n 


2 8-1 
3 6 4 


,and B = 9 -2 
,6 3 


then 


a 1 .b. 1 ={2 8 1)9=2*1 + 8*9 + 1*6 = 80 


r 12 = a 1 .b. 1 =( 2 8 l) —2 =2*7 + 8*(-2) + l*3 = 1 


r 21 =a 2 Jb. i=(3 6 4) 9 =3*1 + 6*9+4*6 = 81 


and 
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ft 2 * ^*2 


' 7" 

4) -2 =3*7 + 6*(-2) +4*3 


Hence, 


For matrix multiplication to be legal, the first matrix must have as many columns as the 
second matrix has rows. This, of course, is the requirement for multiplying a row vector 
by a column vector. The resulting matrix will have as many rows as the first matrix and 
as many columns as the second matrix. Because A has 2 rows and 3 columns while B has 
3 rows and 2 columns, the matrix multiplication may legally proceed and the resulting 
matrix will have 2 rows and 2 columns. 

Because of these requirements, matrix multiplication is usually not commutative. 
That is, usually AB ^ BA . And even if AB is a legal operation, there is no guarantee that 
BA will also be legal. For these reasons, the terms premultiply and postmultiply are often 
encountered in matrix algebra while they are seldom encountered in scalar algebra. 

One special case to be aware of is when a column vector is postmultiplied by a row 
vector. That is, what is 

"- 3 " 

4 (8 2)? 


In this case, one simply follows the rules given above for the multiplication of two 
matrices. Note that the first matrix has one column and the second matrix has one row, 
so the matrix multiplication is legal. The resulting matrix will have as many rows as the 
first matrix (3) and as many columns as the second matrix (2). Hence, the result is 

'-3^ (-24 -6^| 

4 (8 2) = 32 8 

l 7j l 56 U) 

Similarly, multiplication of a matrix times a vector (or a vector times a matrix) will also 
conform to the multiplication of two matrices. For example, 

(8 5Y2^| 

I 6 1 I 8 


9 4 A 5 


is an illegal operation because the number of columns in the first matrix (2) does not 
match the number of rows in the second matrix (3). However, 



Psychology 7291: Multivariate Statistics (Carey) 8/27/98 


Matrix Algebra - 5 


8 5 ^ 
6 1 

9 4 J 


3) 




^8*3 + 5*7) 
6*3 +1*7 
v 9*3 + 4* 7 j 


^ 59 ) 
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and 


(2 7 3) 


5^ 

1 

4y 


(2* 8+ 7* 6 + 3*9 2*5 + 7*l + 3*4)= (85 29) 


The last special case of matrix multiplication involves the identity matrix, I. The 
identity matrix operates as the number 1 does in scalar algebra. That is, any vector or 
matrix multiplied by an identity matrix is simply the original vector or matrix. Hence, al = 
a, IX = X, etc. Note, however, that a scalar multiplied by an identify matrix becomes a 
diagonal matrix with the scalars on the diagonal. That is. 


01 

1, 


01 

4y 


not 4. This should be verified by reviewing the rules for multiplying a scalar and a matrix 
given above. 

Matrix Transpose: 

The transpose of a matrix is denoted by a prime ( A' ) or a superscript t or T (A 4 or 
A t ). The first row of a matrix becomes the first column of the transpose matrix, the 



(2 7 1) 


r 2 8' 

A = 


, and A 4 = 

7 6 


*+ 

00 






vl 4 y 


The transpose of a row vector will be a column vector, and the transpose of a column 
vector will be a row vector. The transpose of a symmetric matrix is simply the original 
matrix. 


Matrix Inverse: 

In scalar algebra, the inverse of a number is that number which, when multiplied 
by the original number, gives a product of 1. Hence, the inverse of x is simple 1/x. or, in 
slightly different notation, x _1 . In matrix algebra, the inverse of a matrix is that matrix 
which, when multiplied by the original matrix, gives an identity matrix. The inverse of a 
matrix is denoted by the superscript "-1". Hence, 

AA 1 = A 1 A = I 

A matrix must be square to have an inverse, but not all square matrices have an 
inverse. In some cases, the inverse does not exist. For covariance and correlation 
matrices, an inverse will always exist, provided that there are more subjects than there are 
variables and that every variable has a variance greater than 0. 
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It is important to know what an inverse is in multivariate statistics, but it is not 
necessary to know how to compute an inverse. 


Determinant of a Matrix: 

The determinant of a matrix is a scalar and is denoted as I A I or det(A). The determinant 
has very important mathematical properties, but it is very difficult to provide a 
substantive definition. For covariance and correlation matrices, the determinant is a 
number that is sometimes used to express the "generalized variance" of the matrix. That 
is, covariance matrices with small determinants denote variables that are redundant or 
highly correlated. Matrices with large determinants denote variables that are independent 
of one another. The determinant has several very important properties for some 
multivariate stats (e.g., change in R 2 in multiple regression can be expressed as a ratio of 
determinants.) Only idiots calculate the determinant of a large matrix by hand. We will 
try to avoid them. 


Trace of a Matrix: 

The trace of a matrix is sometimes, although not always, denoted as tr(A). The trace 
is used only for square matrices and equals the sum of the diagonal elements of the 
matrix. For example. 


( 


tr 
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3 

-1 
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7 

6 
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2 2 

4 

- 5 , 


3 + 6-5 = 4 


Orthogonal Matrices: 

Only square matrices may be orthogonal matrices, although not all square matrices 
are orthogonal matrices. An orthogonal matrix satisfied the equation 

AA* =1 

Thus, the inverse of an orthogonal matrix is simply the transpose of that matrix. 
Orthogonal matrices are very important in factor analysis. Matrices of eigenvectors 
(discussed below) are orthogonal matrices. 

Eigenvalues and Eigenvectors 

The eigenvalues and eigenvectors of a matrix play an important part in 
multivariate analysis. This discussion applies to correlation matrices and covariance 
matrices that (1) have more subjects than variables, (2) have variances > 0.0, and (3) are 
calculated from data having no missing values, and (4) no variable is a perfect linear 
combination of the other variables. Any such covariance matrix C can be mathematically 
decomposed into a product: 
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C = ADA 1 

where A is a square matrix of eigenvectors and D is a diagonal matrix with the 
eigenvalues on the diagonal. If there are n variables, both A and D will be n by n 
matrices. Eigenvalues are also called characteristic roots or latent roots. Eigenvectors 
are sometimes refereed to as characteristic vectors or latent vectors. Each eigenvalue 
has its associated eigenvector. That is, the first eigenvalue in D (d n ) is associated with the 
first column vector in A, the second diagonal element in D (i.e., the second eigenvalue or 
d 22 ) is associated with the second column in A, and so on. Actually, the order of the 
eigenvalues is arbitrary from a mathematical viewpoint. Elowever, if the diagonals of D 
become switched around, then the corresponding columns in A must also be switched 
appropriately. It is customary to order the eigenvalues so that the largest one is in the 
upper left (, d n ) and then they proceed in descending order until the smallest one is in d nn , 
or the extreme lower right. The eigenvectors in A are then ordered accordingly so that 
column 1 in A is associated with the largest eigenvalue and column n is associated with 
the lowest eigenvalue. 

Some important points about eigenvectors and eigenvalues are: 

1) The eigenvectors are scaled so that A is an orthogonal matrix. Thus, A T = A" 1 , and 
AA T = I. Thus, each eigenvector is said to be orthogonal to all the other eigenvectors. 

2) The eigenvalues will all be greater than 0.0, providing that the four conditions 
outlined above for C are true. 

3) For a covariance matrix, the sum of the diagonal elements of the covariance 
matrix equals the sum of the eigenvalues, or in math terms, tr(C) = tr(D). For a correlation 
matrix, all the eigenvalues sum to n, the number of variables. Furthermore, in case you 
have a burning passion to know about it, the determinant of C equals the product of the 
eigenvalues of C. 

4) It is a royal pain to compute eigenvalues and eigenvectors, so don't let me catch 
you doing it. 

5) VERY IMPORTANT : The decomposition of a matrix into its eigenvalues and 
eigenvectors is a mathematical/geometric decomposition. The decomposition literally 
rearranges the dimensions in an n dimensional space (n being the number of variables) in 
such a way that the axis (e.g., North-South, East-West) are all perpendicular. This 
rearrangement may but is not guaranteed to uncover an important psychological construct 
or even to have a psychologically meaningful interpretation. 

6) ALSO VERY IMPORTANT : An eigenvalue tells us the proportion of total 
variability in a matrix associated with its corresponding eigenvector. Consequently, the 
eigenvector that corresponds to the highest eigenvalue tells us the dimension (axis) that 
generates the maximum amount of individual variability in the variables. The next 
eigenvector is a dimension perpendicular to the first that accounts for the second largest 
amount of variability, and so on. 
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Mathematical Digression: eigenvalues and eigenvectors 

An important mathematical formulation is the characteristic equation of a 
square matrix. If C is an n by n covariance matrix, the characteristic equation is 

1C- XII = 0 (1.0) 

where X is a scalar. Solving this equation for X reveals that the equation is a nth degree 
polynomial of X. That is, there are as many Xs as there are variables in the covariance 
matrix. The n Xs that are the roots of this polynomial are the eigenvalues of C. Because C 
is symmetric, all the Xs will be real numbers (i.e., not complex or imaginary numbers), 

although some of the Xs may be equal to or less than 0. The Xs can be solved for in any 
order, but it is customary to order them from largest to smallest. 

To examine what is meant here, let C denote a two by two correlation matrix that 
has the form: 

1 

vP 1, 

Then the quantity C - XI may be written as 


1 p^ 


'X O' 


'l-X p j 

vP ly 


,0 b 


V p 


The determinant is 

|C-Xl|=(l-X) 2 -p 2 

So the equation that requires solution is 

(l-X) 2 -p 2 = 0. 

This is a quadratic in X. If we had three variables, it would be a cubic; and if there were 
four variables, it would be a quartic, etc. Solving for the quadratic gives 

X = 1 ±p 

The largest root depends on the sign of p. For p > 0, then X 1 = 1 + p and X 2 = 1 - p. 

For each X, one can define a nonzero vector a, such that 

(C-Xl)a = 0. 

The 0 to the right of the equal sign denotes a vector filled with Os. Any number in a that 
satisfies this equation is called a latent vector or eigenvector of matrix C. Each eigenvector 

is associated with its own X. Note that the solution to a is not unique because if a is 
multiplied by any scalar, the above equation still holds. Thus, there are an infinite set of 
values for a, although each solution will be a scalar multiple of any other solution. It is 
customary to normalize the values of a by imposing the constraint that a' a = 1. A latent 
vector subject to this constraint is called a normalized latent vector. 

Taking the two by two correlation matrix with p > 0, then 
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(C-Al)a = 


'1-A P) 

( a 0_' 

V p 

UJi 


(l-?i)a 1 +p a 2 A 
(1- A,) a 2 +paj 



or, by carrying out the multiplication, we find 

(l-A)^ +p a 2 =0 
(l-A )a 2 +p a 1 =0 

Now take the largest eigenvalue, 1 = 1 + r, and substitute. This gives 

p (a 2 -a 1 ) = 0 
p (a 1 -a 2 ) = 0 

Thus, all we know is that a 1 = a 2 . If we let a 1 = 10, then a 2 = 10; and if we let a 1 = -.023, then 
a 2 = -.023. This is what was meant above when it was said that there were an infinite 
number of solutions where any single solution is a scalar multiple of any other solution. 
By requiring that a'a = 1, we can settle on values of a 1 and a 2 . That is, if a t 2 + a 2 = 1 and a 1 = 

a 0 = a, say, then 2 a 2 = 1 and a = J^5. So the first eigenvector will be a 2 by 1 column vector 
with both elements equaling Js . 

For the second eigenvector, we substitute 1 - p for A. This gives, 

p(fl! + a 2 ) = 0 

Consequently, a 1 = -a 2 . One of the a's must equal -J3 and the other must equal -J^5. It is 
immaterial which is positive and which is negative, but it is a frequent convention to 
make the first one (flj) positive and the second negative. 

Note that the normalized eigenvectors of a two by two correlation matrix will 
always take on these values. The actual value of the correlation coefficient is irrelevant as 
long as it exceeds 0. 

Can you figure out the eigenvectors for a two by two matrix where p < 0? 



